Quantum Hall effect in honeycomb lattice under strain 
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We investigate the possibility of the quantum Hall effect in the deformed honeycomb lattice. 
Static strain in the honeycomb lattice preserves the time reversal symmetry (TRS), and can be 
viewed as a valley dependent gauge field in the continuum limit. We calculate the response of the 
system to the valley magnetic field induced by strain and compute the valley Hall conductance. It 
is shown that the system exhibits a nontrivial response (nonzero Hall conductance) even in time 
reversal invariant (TRI) samples. We argue that a quantized current will be induced if we subject 
the system to a valley electric field by applying time dependent strain. A theory for the low energy 
excitations of the systems is developed. We show that the mutual Chern-Simons action describes 
the low energy theory of this system and also discuss the necessary details for observing such a 
behavior with a proposal for experiment. 
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Introduction. — The emergence of the fractional charge 
and statistics in quantum Hall systems, has spurred a 
broad interest in topological states of matter. In semi- 
conductor quantum Hall liquids, strong magnetic field at 
integer and fractional filling fractions can generate quan- 
tized off-diagonal conductivity [T]. The experimental 
discovery of the quantum Hall effect (QHE) in mono- 
layer graphene, has attracted a lot of attention to the 
properties of two-dimensional (2D) Dirac fermions in the 
presence of magnetic field [2] . Haldane in a seminal pa- 
per [3] showed that in order to have integer quantum Hall 
effect on honeycomb lattice, the presence of the external 
magnetic field is not necessary. He suggested a model 
in which the time reversal symmetry (TRS) breaks be- 
cause of the pattern of the hopping integrals, and the Hall 
conductance does not vanish despite the vanishing total 
magnetic field. Accordingly, it is believed that the time- 
reversal symmetry breaking by itself can potentially lead 
to integer QHE. This has been recently verified for par- 
tial filling fractions as well j3]. It is also well established 
that topological insulators which preserve TRS exhibit 
quantum spin Hall effect [3. In this paper we dis- 
cuss that static strain which does not break time reversal 
symmetry induces a nonzero valley magnetic field whose 
direction changes as we travel from one valley to anther 
in the honeycomb lattice. We show that this valley mag- 
netic field induces carrier density with a quantized valley 
Hall conductance. We discuss through the Laughlin ar- 
gument and the Chern-Simons theory that a time depen- 
dent strain which induces a nonzero valley electric field as 
well, generates a nonzero current in the system. It is also 
shown that the low-energy behavior of this system can 
be described by the mutual U{1) x U{1) Chern-Simons 
action j^. 

In the following, we discuss a possible realization for 
suggestion for charge quantum Hall effect (CQHE). In 
the honeycomb lattice with massive Dirac fermions, for 
each spin degree of freedom, the low energy description is 
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FIG. 1: (Color online) (a): Upper hemisphere represents the 
Berry curvature for Hk,s which leads to Ck = Ck,^ + C'x 4 = 
1. Similarly, lower hemisphere represents if 'a. According to 
Eq. (1), the pseudo-magnetic field induces the same charge 
density at K and K' valleys, (b): The sign of the spatial 
current along the a:-axes is determined by CEy. Hence, K 
and K' contribute the same current. 

given by the sum of two independent sub-Hamiltonians. 
We can tune the structure of these sub-Hamiltonians to 
have opposite Chern numbers. On the other hand, it 
has been discussed in the literature that strain on hon- 
eycomb lattice can induce an effect similar to that of the 
magnetic field, which is called pseudo- magnetic field 
Basically, strains changes the bond length between atoms 
and due to the particular band structure and lattice sym- 
metry of the honeycomb lattice, strain creates pseudo- 
magnetic fields acting on the charge carriers. Because 
strain does not break TRS, the two sub-Hamiltonians 
feci opposite pseudo-magnetic fields. Massive pseudo- 
magnetic fields raised by highly strained nano-bubbles 
that form when graphene is grown on a platinum sur- 
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face have been recently measured by STM spectroscopic 
measurements from the spectral peak spacing [9]. The 
spectral peak indicates Landau levels created by the local 
pseudo-magnetic fields. Accordingly, strain engineering 
can prepare the desired system. 

Since spatial current of edge states in the QHE is de- 
termined by the product of the Chern number and the 
sign of the pseudo-electric field [TD], both bands con- 
tribute the same current and therefore the total cur- 
rent is nonzero. More specifically, let Hi and H2 be 
Hamiltonians with Chern numbers Ci and C2, respec- 
tively. We assume that Hi and H2 are subject to Bi 
and B2 pseudo-magnetic fields. Then the conductiv- 
ity is determined by Ji^x — \GiEi^y where i G {1,2}, 
and Ji^x is the x— component of the current induced by 
the pseudo-electric field Ei^y. The current thus depends 
on CiEi^y. Requiring the time reversal invariance of 
Hi (B H2, we obtain C2 — — Ci and B2 — — Bi , whence 
A2 = — Ai. Therefore E2,y ~ —Eiy, and the total cur- 
rent Jx = Ji.x + J2,x and the quasiparticle density are 
given by: 
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J. = \{CiEi^y + C2E2,y), n = \{CiBi,z + C2B2M 

It means that those circumstances are sufficient to have 
a nonzero total edge current. There is also another way to 
justify the above claim. Within the continuum model. Hi 
and H2 are independent and are invariant under global 
f/(l) gauge transformations. The corresponding Noether 
currents are Ji and J2. Therefore, both J+ = Ji + J2 
and J_ = Ji — J2 are conserved quantities jlT]. J+ is 
the electromagnetic response of the system, while J_ is 
the response to the pseudo-electromagnetic field. To this 
end, note that everywhere q^A^ appears in the Hamilto- 
nian where qi is a particle charge. We can safely assume 
A2 = Ai and 92 = —q_i- Consequently, A, acts like ordi- 
nary electromagnetic gauge field which couples to oppo- 
site charges at two valleys. This is a way to detect the 
conserved current J_, which can in principle be nonzero. 

In Ref. [12] , a condition to determine a zero-field QHE 
in massless graphene and quantum pseudo-spin Hall ef- 
fect like topological insulators in the presence of pseudo- 
spin-orhii coupling have been proposed. Nonetheless, 
since the two generated edge states here circulate in op- 
posite direction, the net current would be zero; thus in- 
stead of a CQHE a pseudo-spin QHE would emerge. 

Physical realization. — Let us consider the Hamiltonian 
of whether undoped graphene or an artificial 2D honey- 
comb lattices in the continuum approximation [31 113j . In 
that limit it is given by the sum of the following two 
independent massless Dirac Hamiltonians {h — 1): 



where Q = {K, K'} and a^'V-^ are Pauli matrices, K 
and K' refers to two dirac points and s is the spin degree 
of freedom. The sign +(— ) refers to K{K'), respectively. 
Because we are interested in the Hall conductance of the 
system, we need to compute the Chern number of the 
ground state. We require the resulting Hamiltonian i.e. 
Hk ® Hji' to be time reversal as well as spin rotation 
invariant. To this end, we do need an opening gap mech- 
anism in the system. In this case, the mass term for 
both valleys should be equal and independent of the spin 
of electrons. Therefore, we should have: 
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It is easy to check that the above Hamiltonian neither 
breaks the time reversal nor the spin rotation symmetry. 
A generic noninteracting two band Hamiltonian can be 
written a.s H ^ -EfcV''^(^) (doik)! + d{k).a^ ip{k). For 
every fc, the term in parenthesis resembles the Hamilto- 
nian for a spin 1/2 in a magnetic field given by dk, and 
therefore, ip{k) aligns itself with the direction of dfc. Thus 
the Hamiltonian maps the k space (2D torus) to a unit 
2D sphere spanned by dk- The Chern number can be 
viewed as the winding number of this mapping [3 110) . 
It can also be computed through C — where 6*5 is 

the Berry phase of the ground state along the path 

given by (0,0) ^ (27r, 0) ^ (27r,27r) ^ (0, 27r) (0,0). 
Using this recipe, it can be shown that the Chern number 
of the massive Dirac Hamiltonian of the following form 

H ^ J d^xi;^K,a i^<^'' + iVi^^'dx + imcr^dy) tpK,a (4) 

is given by C = = ^sgn (771772 m) |[7j. Therefore, we 
have: 
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C'k' — ^ Ck',ct 



(5) 



Hq^s = VF / d'x^P^Q ^{x) (±za"c), + la^dy) i^q^x) (2) 



Because we require the system to be time reversal and 
spin rotation invariant, the mass term should be inde- 
pendent of the valley degree of freedom as well as the 
spin of quasiparticles. One possible way to generate such 
a mass term is to introduce staggered chemical potential 
i.e. on sublattice A we have fiA = m and on sublattice B 
have /is — —m. This pattern of the chemical potential 
generates a mass term as we desire. This mechanism, for 
instance can be realized experimentally in graphene by 
chemical doping [TH [TS] . 

Another mechanism to generate the mass term is to 
take strong on-site repulsions, U, into account e.g. by 
adding the Hubbard term. The system is no longer 
graphene in that case, but still lives on honeycomb lat- 
tice |16[ . In the rest of this paper we will demonstrate 
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that our result depends only on the honeycomb lattice 
structure and applies to this case as well. Several au- 
thors including one of us, have studied the properties 
of the Hubbard model near the Mott transition. Both 
numerical and theoretical studies have reported the exis- 
tence of gapped spin liquids p!7H^ . In this state of mat- 
ter, the fermionic degree of freedom can be modeled by 
two gapped Dirac Hamiltonians for each spin with a mass 
term of the form we used in the previous section. There- 
fore, slightly above the Mott transition {U > Uc), but 
still below the onset of anti-ferromagnetic order, a mass 
term will be added to the system which respects all sym- 
metries of the lattice including time reversal, spin rota- 
tion, parity, 60*^ rotation and translation symmetries |21j . 

Strain engineering and induced pseudo-magnetic 
field. — It is interesting to investigate the effect of strain 
on defect-free honeycomb lattices due to smooth elas- 
tic deformations. Strains can create induce gauge field 
and of at least 10% can be induced in the system with- 
out damaging its structure [22] . Importantly, due to the 
gauge field of strains, quasiparticles near K point feel 
a pseudo-magnetic field that is opposite to the pseudo- 
magnetic field that quasiparticles near K' feel, in agree- 
ment with the requirement of the time reversal invari- 
ance. Therefor we replace Hk and Hk' with the follow- 
ing Hamiltonians, respectively (e — 1): 



HK,a — 



Hk' 



^cjy (idy-Ay) 



v^a^ {idy 



ll^K'M (6) 



A connects to the strain of the sample like 
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where 



is the strain distribution, 



d\og{t) / d\og{a) w 2, a w 1.4 A is the lattice con- 
stant, t ~ 2.1 eV is the nearest neighbor coupling energy, 
and c = ^/?,^^/{2n) is a constant and n « 0.56 [SF. $0 
is the quantum unit of magnetic flux. We will propose 
a feasible way to generate proper strains and pseudo- 
magnetic fields on strained graphene nanoribbons. 

Now at this stage, we use auxiliary gauge field, oq.^ 
to write the current of quasiparticles living on the 
Dirac cone at valley Q (where Q = {K,K'}) as Jq = 
^(-'^'"''dyaQ^p {h — 1, e = 1). The reason is that 
d^jq = is naturally satisfied using that definition. We 
can obtain the effective action for the gauge field degree 
of freedom, which describes the low energy theory of the 
system under consideration, by integrating out fermions. 
Because the Chern number for both Hk and Hk' is non- 
zero, it is known that the Chern-Simons action describes 



the effective low energy theory for both, and therefore 
the Lagrangian corresponding to Hq is: 

Cq = ±^e^-PaQ,^d,aQ^, + ^^e^^^" AQ^.d^aQ^, (8) 

where ± denotes the valley indices and = 
[Aq, Ax, Ay). Without any loss of generality, we assume 
Aq = Q for simplicity. The second term in Eq. (jsj which 
is JqAq^p, is the minimal coupling between the induced 
gauge field and the current of quasiparticles. The La- 
grangian of the whole system is then C ~ Ck + Ck'- 
The corresponding action, S = J d^xdt C, is known as 
mutual J7(l) X J7(l) Chern-Simons action. Although both 
Ck and Ck' break the TRS, their sum does not. This 
statement becomes obvious by considering the following 
transformations under the TRS: 
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Ak.p 
{a-K.o, ^k) 



K' 

{-aK'fli 



SLK' 
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The Lagrangian in Eq. ([8]) leads to the following equa- 
tion of motion: 

fg - ^^""d.aQ^p = ±l-e^-Pd.AQ^p (10) 

restoring h and e, the total current is 



2e2 
ti 



(11) 



finally, the pseudo-magnetic and the pseudo-electric fields 
due to strain lead to the following induced currents 



y,K, 



Jy = -j^Ex^K, 



2e' ^ 



Bk (12) 



It is worth mentioning that the induced current cre- 
ates a real vector potential, Aem from Maxwell equa- 
tion V X V X Aem = MoJ- This vector potential couples 
to the current of quasiparticles in the bulk j, through 
— e J d^x Aem ■ j term in the Lagrangian. Since both 
Aem and j change sign under TRS, then the excess term 
is TRI and so is the total Lagrangian. It is worthwhile 
mentioning that Aem is weak and since Ck + Ck' = 0, 
it does not generate edge current. Consequently, it can 
be neglected up to the first order of approximation. 

One important consequence of these results is that 
the charge quantum Hall effect due to strains can serve 
as a new probe to detect topological insulator phase 
in graphene. In topological insulators, the charge cur- 
rent is zero, even in the presence of pseudo-electric field 
induced by strains. To see this, consider the Chern 
number of spin up electrons at valley K. We have 
Ck.^ — —Ck.-I — \- Therefore the current at K, is 

Jk = Jk,-\ + JK,i = X ^'^K,-\ + Ck,x) Ek = 0. Similarly 
we have Jk' = and therefore J = Jk + Jk' = 0. Con- 
sequently, we can measure the induced charge current in 
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a time reversal insulating graphene sample in response 
to strains. For topological insulators, on the other hand, 
this current will be zero, while for trivial insulators it is 
nonzero. 




FIG. 2: (Color online) (Left): Proposed strain for a sample 
of graphene, lOOA x WOA with Pb = 0.002A"^ The sample 
size is just used for demonstration. (Right): Direction of the 
charge current for different Dirac points in different edges. 

We propose an experimental setup to measure 
CQHE based on recent developed techniques to grow 
graphene on a thin cuprates showing outstanding flex- 
ibility [53] subject to chemical doping using hydrogen- 
passivated [2]. The atomic displacements are given by 
iu^,Uy) = {-2Pb xy + PE{t)y^,PB{x^ + 2/^)) in the 
Cartesian coordinates (Fig. [2| where Pb and Pe (t) are 
geometric parameters with unit of 1/length and Psit) 
can be a time-varying periodic function, for instance 
PE{t) = PBCOs{ojt). The correspondence strain leads to 
a uniform magnetic field = ^(pQcP Pb / a and pseudo- 
electric field Ey = —2(j)Qc(3uJsm{ujt)y/a. The induced 
scalar potential, which is proportional to u^x + ""j/y, 
vanishes due the special introduced strain. Consider a 
strain geometry Pg = 0.5 {^m)~^ in a sample with 
size 0.4 /im x 0.4/im, the pseudo-magnetic field yields 
Bz « 9.0 T which is large enough to observe the quan- 
tized quantum Hall Landau levels. Due to the finite 
width, there are finite edge modes in the zigzag bound- 
ary |27| . The zigzag boundary condition used here do not 
mix the K and K' valleys, leading to a clear distinction 
of the edge modes. Besides the edge modes, an exter- 
nal electric field can not induce pure current along the 
horizontal boundaries, however the pseudo-electric filed 
constructed by —dA{t) /dt creates easily a voltage along 
the boundaries demonstrated in Fig. [2] Importantly, 
using the parameters and uj — 60 AIHz, the maximum 
pseudo-electric field reaches 50 V/m and then would 
be about 4 x 10"'^ A/m. This induced current can be 
measured by experiments. 

In summary, we have shown that CQHE can be seen 
even without the TRS breaking in a strained honeycomb 
structures which can be either artificial honeycomb lat- 



tices or graphene nanoribbons. We have introduced stag- 
gered chemical potential to open a gap in the band struc- 
ture together with a varying induced gauge field created 
by strain on the system. Furthermore we have shown 
that, even though that the whole system does not break 
the time reversal, the total current in the edge of the 
system is nonzero and we then observe CQHE. Finally, 
our theory can be investigated by an experimental setup 
that we have proposed. 
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